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Abstract. Topological aspects of graphene are reviewed focusing on the massless Dirac 
fermions with/without magnetic field. Doubled Dirac cones of graphene are topologically 
protected by the chiral symmetry. The quantum Hall effect of the graphene is described by 
the Berry connection of a manybody state by the filled Landau levels which naturally possesses 
non-Abelian gauge structures. A generic principle of the topologically non trivial states as the 
bulk-edge correspondence is applied for graphene with/without magnetic field and explain some 
of the characteristic boundary phenomena of graphene. Q 



1. Introduction 

Graphene as a playground of massless Dirac fermions has a long history of theoretical 
study [U El [3]. However its experimental realization is a real surprise since the massless Dirac 
fermion itself is topologically non trivial especially under a magnetic field [H [5]. It implies that 
we have a chance to observe exotic topological phenomena in a real world, which have been 
confirmed experimentally for the past few years. 

The quantized Hall conductance is one of the most well-known topological objects as the 
Chern numbers of the Berry connection which describes the bulk (without boundaries). It 
is also described by the edge states, supplemented by the Laughlin argument, which gives 
a different topological number. These two topological quantities are closely related as the 
bulk-edge correspondence [6j. This bulk-edge relation has been realized as a universal feature 
of topological states in several different quantum systems such as superconductors, quantum 
(spin) Hall effects O [S], gapped quantum magnets and photonic crystals[10j. The edge states 
in a magnetic field are chiral fermions in a sense, they have specific direction, and can not 
be destroyed as far as the bulk is gapped. It is the topological stability of the edge states. 
Graphene under a magnetic field is one of the typical topological states where the bulk-edge 
correspondence of the Dirac fermions plays a fundamental role. Localized zero modes near the 
zigzag boundaries (Fujita states) [TT], which have been measured by the STM experiments [T2] . 
and their generalization under a magnetic field, are also governed by this topological principle [6]. 
Further appearance of two Dirac fermions with opposite chirality at the K and K' points are 
not accidental but is topologically protected by the chiral symmetry which is a two dimensional 
analogue of the Nielsen-Ninomiya theorem in four dimensions [13J. 

Also the Landau level of graphene is quite special reflecting its singular dispersion as of the 
Dirac fermions. Then the quantum Hall effects of it is anomalous |141 \TE[ 116]. Further Landau 
degeneracy of the characteristic zero energy Landau level are topologically protected against for 
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some class of randomness by the index theorem. It is fundamental for free-standing single layer 
graphene where its intrinsic disorder is a gauge field fluctuation [TTJ [18]. 

Based on our works, these topological phenomena in graphene physics will be explained 
intuitively. 



2. Topological stability of the Dirac cones and the n = Landau level 

2.1. Chiral symmetry and doubled Dirac cones 

Let us consider a single orbital tight-binding hamiltonian of graphene with magnetic field |16j 
(J = (ill 32) is a t wo dimensional coordinate with the unit translations ei i2 and (j) is a total flux 
per hexagon in flux quantum $0 = h/e, See Fig(T|) 
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where we give two expressions using two different unit cells that are compatible for zigzag 
and bearded edges. When <p = 0, it is written in the momentum representation as H = 
J^(k)h(k)c(k), 
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where c^(fe) = (cj (fc), cl(fc)), c a (j) = J e lk ^c a {k). The energy dispersion is given as 
e(fc) = ±|A(fc)| and vanishing momenta kp, A(krj) = give massless Dirac cones if exist. 




Figure 1. The unit cell of the 

graphene. The red ellipse is consistent Figure 2. Trajectories of A(fc) parametrized by 

with the zigzag edges and the blue one fci for various values of k 2 . (a) for A = A z and 

is for the bearded one. Two primitive (b) A = A B (i = — 1). Curves that encloses the 

translation vectors are also shown. origin is drawn in red. 



Since the Brillouin zone, T 2 = {(kx, kz)\ki, G [0, 2ir\), is a two-torus by identifying 
k\ = ^ ki = 27T and A; 2 = fc 2 = 2ir, it does not have boundaries. Then considering 
ki as a parameter, A(fe) makes a closed loop in the complex plane for each fixed /c 2 . This loop is 
modified by changing fe 2 from to 2ir and comes back to the original one (FigJ2J). When the loop 
cuts the origin, it gives the Dirac cone. Then once the loop cuts the origin, the situation is stable 
against for small but finite perturbation. This is the topological stability of the Dirac cones. 
Here we need the 2x2 hamiltonian is characterized by single complex parameter A ( diagonal 



parts are zero). This is the chiral symmetry, that is, there exists {h, 3 ~f} = 0, -y 2 = a®: 2x2 
unit matrix. (In this case, 7 = 03. ) Also the doubling of the Dirac cones is again clear since 
the origin cutting the loop from the outside to the inside inevitably followed by cutting from 
inside to the out (See Figj2] and FigU]). This simple observation guarantees that the number 
of Dirac cones are always even. Note here that the chiral symmetry is related to the bipartite 
structure of the honeycomb lattice, that is, the hoppings are only between the two sub-lattices 
o and • 

On the other hand, as for generic semiconductors including graphene, let us assume the 
effective hamiltonian of the valence and the conduction bands, h, is chiral symmetric, that is, 
it does anti-commutes with some chiral matrix 7, {^1,7} = 0. It implies h is traceless and thus 
one can expand it by the Pauli matrices as h(k) = R(k) ■ cr with three real parameters R±(k), 

( h ' i k > \ 

R2(k) and Rs(k) that makes a vector R(k) = Rz(k) . The momentum dependent energy 

V Rs(k) J 

gap is E g {k) = 2|i2(fe)|. Then the zero gap condition is R = 0. Writing the chiral matrix as, 
7 = n 7 • cr (n 2 = 1), the chiral symmetric condition {h, 7} = (R • n 7 )o"o + i(R x n 7 ) • cr + (Rt± 
n 7 ) = 2(R-n^)ao = 0, implies i?-n 7 = 0, that is, R is always in the plane M 2 (n 7 ) whose normal 
vector is n 7 . Since the Brilluine zone T 2 is a 2-dimensional closed surface, the image R(T ) 
is also a closed surface (like a balloon) in 3-dimensions. Then the condition R(T 2 ) C M 2 (n 7 ) 
means the closed surface is collapsed on the plane (a rubber balloon without air on the desk). 
When the collapsed image includes the origin, it gives the Dirac cones. Then the topological 
stability of the Dirac cones and the doubling of them are clear |1 6\ \19\ I20j . (See Figa3J) It is 
the 2-dimensional analogue of the Nielsen-Ninomiya theorem in the 4-dimensional lattice gauge 
theory |13j . A relation between the 4D graphene and chiral fermions is recently discussed as 
well [21]. 




Figure 3. Topological stabilities of the Dirac cones. The image of the Brilloine zone R{T 2 ) is generically 
a surface in three dimensional space R, which is collapsed into the plane which is normal to the n 7 , when 
the hamiltonian is chiral symmetric as {H, 7} = 0, 7 = cr • n 1 . 



2.2. Chiral symmetry and the n = Landau level 

Now let us expand the hamiltonian h(k) with chiral symmetry ({h, 7}) near one of the doubled 
zero gap momentum, say, k^. Then the effective hamiltonian in the lowest order in 5k = k—ko is 
written as h ~ fac" = Sk\<T-X+5k,2<r-Y where X = d^R and V = dt 2 R which are perpendicular 
to the n 7 . It implies the chiral symmetry of the effective hamiltonian, {hc,j} = 0. One may 
define the chirality of the effective hamiltonian xd = if as n-y = XdX x Y j{cK) 2 where the 
effective light velocity (c > 0) is defined as c 2 = \X x Y\/h 2 . When the vector, (X ,Y ,n 7 ), 
forms a right handed triple, the chirality, xd is +1 and, if left-handed, xd = — 1- Then the 
chiral matrix is written as 7 = XD& • (X x Y)/(ch) 2 . 

By the inverse procedure of the quantization, H5k — > p — eA = 7r, p a = —ihd a , (a = x, y), we 
have a real space form of the effective hamiltonian with a magnetic field, B = d x A y — d y A x > 0, 
[tt x , TT y ] = iheB as 

h c = ^[(ff-lK + ^TjirJ 

The spectrum and the wave functions are determined by considering its square as 

h 2 c = h- 2 [X 2 n 2 x + Y\ 2 + (X-Y)(7T x Tr y + Vx ) + i(X xY)-a[TT x ,TT y ]] 
= c 2 (7r^H 7t)<to — XD(eBhc 2 )j 



where 
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By identifying as i = c 2 , the first term of he is considered as a hamiltonian of the standard 
parabolic electrons with anisotropic masses (See |Appendix A ). Here the cyclotron frequency cue 



is identified as ojq = eB/m = 2eBc . Also note that the chiral matrix, 7, commutes with the 
squared hamiltonian h 2 c . 

Then writing a normalized standard Landau state with energy, tvjj c {n + 1/2) = 2heBc 2 (n + 
1/2), n = 0,1, 2, • • • , as tp n (x,y), the eigen state of h 2 c is written as = ^>±ip n (x,y) 
where ip± is an eigen state of 7 with eigen values ±1 respectively. It satisfies h 2 -,^^ = 
{eBhc 2 ){2n + 1 =F xd)^u- ^ implies that the effective hamiltonian he always has a zero energy 
Landau level as 



v XD ( x >y) = v XD Mx,y) ■ h c v° XD {x,y) = 

which has the chirality xd- 

As for the non zero energy Landau levels, h c is doubly degenerated at e 2 , e n = cV^neBh, 
spanned by ^t, XD = f+XD^nix^y) and *~^f = (p- XD ipn-i(x,y), which diagonalize he without 
degeneracy with the energies ±e n . 



2.3. Chiral symmetry and Fermion doubling 

Up to this point, we have discussed the effective theory of the chiral symmetric zero gap 
semiconductor near one of the gapless momentum. 



As discussed the chiral symmetry of the global 
hamiltonian h, {h(k),~f} = 0, (h = R ■ cr ) requires 
that the three dimensional vector R(^k) is always 
perpendicular to the n 7 and on the plane R 2 (n 7 ) cutting 
the origin R = ( R(k) G M 2 (n 7 )). Then one can 
repeat the discussion for A(fc) in Figj2]for the contour RQ 
on the plane M 2 (n 7 ), that is, considering a closed curve 
characterized by A>2, C(k2) = {k\k\ £ [0,2-7r]}. This curve 
is on the plane M 2 (n 7 ) and comes back to the original one 
by changing k\ form to 2ir. When it cuts the origin, it 
gives the Dirac cones. The doubling of the Dirac cones is 
also clear as before (Figs 12] and HJ. Also as for the chirality 
of the Dirac cones, xd, the paired Dirac cones, D and D' , 
have a reversed chirality as 

= 




Xd + Xd> 

which is also clear from the FigJU Then the chirality of the 
n = Landau level is also reversed. As for the graphene, 
the Dirac fermions at K and K' do have the reversed chirality. 



Figure 4. R(k) on the plane K 2 (n 7 ). 
(a) When the curve cuts the origin 
which defines the Dirac cone. (b) 
Another Dirac cone which is paired 
with the one in (a). The chirality is 
reversed. 



2.4- Aharonov-Casher argument 

As discussed the n = Landau level is an eigen state of the chiral operator for each Dirac fermion 
which is a characteristic feature of the chiral symmetric massless Dirac fermions on lattice. It 
is realized in graphene. This complete degeneracy of the Landau level is not only for the pure 
system and it does persist for some class on the randomness [?]. There is well-known discussion 
by Aharonov-Casher [22J which is considered as a direct demonstration of the index theorem 
where the degeneracy is only determined by the total flux passing through the system as far as 
the effective description by the Dirac fermions is allowed. Here the magnetic field B may not 
be uniform. We compactly describe it extending the discussion to our anisotropic case. Since 
the anisotropic mass matrix 5 discussed before is real symmetric and its determinant unity, it is 

V, where 

= iheB], 



diagonalized by the orthogonal matrix V as (See also appendix), H = 

£,x,£,Y > 0, £x£y = 1- Then defining canonical momenta II = Vtt, [ILx, Ily] = [■ 
h 2 c is factorized as 
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Then writing the vector potential as Ax = —£,y9\ 
and Y = Y/VZy ( Vlxdx = d 
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and defining X = X/y/£ x 
<p XD if) satisfies, 



tpy ■ Then writing as ipo = e 
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/, the function 



/ satisfies {d x + id Y )f = 0. This is written by defining z = X + iY as d z f = 0. Thus the 
function / is an entire function of z in the whole complex plane z G C, that is, polynomials. Also 
note that the function 4> needs to satisfy B = dxAy — dyAx = (£,xd x +£,Yd Y )4> = (d x +d^.)(j>. It 
implies <j)(X, Y) = J dXdYG(X-X', Y-Y')B(X', Y') where G(X, Y) = ± log r 2 = X 2 +Y 2 
and ro is a constant. When assuming the magnetic field is only nonzero in the finite region, 
we have for the asymptotic behavior in the limit r — > oo as <f> r — ? J^l°g(^) where 
$ = J dXdY B = j dXdY B is a total flux. Then we have for the asymptotic behavior of 



the zero mode as ip r —? f( z ^JL^ <t>o , it implies that the number of the degeneracy is 



3. Hall conductance and Berry connection of Dirac sea 

Geometrical phase of the quantum mechanical states are of fundamental interest for physical 
society today. The quantum Hall effect and its time reversal invariant analogue as the quantum 
spin Hall effect are the typical examples. 

Now let us start from the Niu-Thouless-Wu formula for the Hall conductance [23J. It reads 
that the Hall conductance of the manybody state is given by the average over the twisted 
boundary condition as 
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where 8,, 



d%i (a* 



x,y) and where \G) = \G((j) , 4> 2 )) is a unique gapped ground state of the 

hamiltonian, H(4> 1 , <f> 2 ), with a twisted boundary condition specified by e^ 1 and e 1 ^ 2 for each of 
the two translational directions as H\G) = E\G). 

This C is the topological number (Chern number) and is intrinsically integer as explained 
below. This topological property is clearly demonstrated by defining the Berry connection ( a 
formal form of the vector potential) as 



A 



{G\dG) = {G\dM)d^ = d<f> ■ A, 



A = (G\VG) 



where <f) = (0 1 ,0 2 ,O) and V = (81,82,83). Note that, since the phase of the state \G) is 
completely arbitrary, one may take a different choice as \G') = \G)uj, w = e' 9 , £ K, which gives 
a different connection A' 

A 



(G'\dG') = 
A' + u' 1 



dcf) ■ A' . They are related with each other as 
A' + id6, A = A' + iV9 



dco 



This is the gauge transformation similar to the case of the electromagnetic field governed by the 
Maxwell equations. (As for the differential form, see [25].) This is called Abelian since uj is just 
a complex number \uj\ = 1. 

Since the twisted boundary condition of e l2lT is equivalent to the periodic boundary condition 
without any twist, the integral region T 2 = {(<f> , (p 2 )\(j)^ £ [0,27r]} is identified as the two-torus. 
Then the Chern number of the manybody state is written as 



C= — I dS ■ rot A = — / dA=— I 7 = 
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where dS = 
connection A is 
T' = dA' = dA since d(oj~ 1 dui) 



i(j) 2 (0,0, 1) and the field strength is defined as T = dA = (dG\dG). The 
gauge dependent but the Abelian field strength T is gauge invariant as 
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dujuj' 



-(uj^duj) 2 = (See|Appendix B[). 



To understand the quantization of the Hall conductance, 
we need to specify the gauge [26J. Here we discuss this gauge 
fixing generically|27j. Since the two-torus T 2 is boundaryless, 
the Chern number is always vanishing if one is allowed to use 
the Stokes theorem globally C = j T2 dA = Jg T 2 = A = 0. To 
have a nonzero Chern number, the single global gauge is not 
allowed and one needs to use two different gauges, at least, 
say, Ai and A2- Let us first mention that the state |G) 
is gauge dependent but the projection into the ground state 
space P = \G)(G\ = \G'){G'\ is gauge invariant. Then taking 
arbitrary (fixed) state |Ti), a normalized and gauge fixed state 
Id) is given as |G X ) = P|ri)/V^t where iVi = (Ti|P|Ti) = 




Figure 5. Singularity of the 
gauge and the two patches to 
avoid it. 



|(G|Ti)| 2 > is a gauge invariant normalization constant. It is 
clear that this gauge becomes singular when the normalization 

Ni is vanishing at some point in T 2 (N% = 0). Near this singular point 4>\, one needs to use the 
other gauge, say IG2) using another arbitrary state IT2) that is generically regular in the region 
R ± that includes fa, N 2 ((f>) ^ 0, (j> e R\. 

The gauge transformation between the two, A\ = Ai + 1621 is given as \G\) = \G 2 )oj 2 i, 
W21 = (T2\P\T\) / N1N2 = e 10 ' 21 . When these two gauges are enough to span the whole <fi space, 
the Chern number is written as 

C = ^-( [ dA 2 + I dA^] = [ (Ai -A 2 ) = ^- I d9 2 i 



27ri V JRi ' Jt 2 \R! ) 27ri JdRi 27T J OR! 

Here the reason the Chern number is integer is clear since the gauge transformation oj 2 i is single 
valued on the boundary dR\. 

When the system is non-interacting, V = 0, one may go further. In this case, the ground 
state of the many body system, (M particle state), \G) is constructed from one particle states. 
It is explicitly expressed by writing the hamiltonian as H = c^(h)c where c' = (c\, • ■ ■ , cn) and 
N is a total number of system sites. The ground state is given by the filling the one particle 
states below the fermi energy ep as (em < £f) 

M 

\G) = |*) = n ct ^l°>> hipt = e e ij> e , e e < e t ; (£<£'), $\^v = hv 
£=1 

where the many body state is labeled by the filled one particle state, that form a multiplet, 
\l> = (-01, • • • , i/>m)!2E1 ESI ED]- The Berry connection of this non- interacting ground state has 
an interesting form by the non-Abelian Berry connection Am which is defined as 



A = {g\ V(ct0i)---(c^)----(cV M )|O) = Vdet*t(^i,--- ,d^ e ,--- , if>M ) = Tr A 

* • * * A/7 
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Y M = -- 
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Here the choice of the one particle states below the fermi energy is arbitrary and one may further 
allow to mix up among them as 

* = *'cj, u^u = E M 

which results in the same many body state = \^)' detcj since u € U{M) and detcj £ U(l). 
This Am defines a non-Abelian Berry connection. Its property under the gauge transformation 
is given as[3TJ 127] 

Am = A' M u + uj -1 dw 

This is the gauge transformation of the non Abelian matrix valued connection spanned by the 
filled one particle states. Correspondingly the field strength (matrix) is defined as 

F M = dA M + A 2 M = Y^ F^d^df, = 3^ A v - d v A^ + [A M , A v \ 

fJ,<V 



It transforms as Fm = u 1 F' m uj. This non Abelian field strength is directly related to the field 
strength of the manybody state 

T = Tr Fm 

which gives the TKNN formula when each Landau levels are well separated |32j. It transforms as 
T = Tr Fm = J 7 ' = Tr F' M , which is consistent with the physical observable a xy is independent 
of the choice of the one particle state. 

The advantage of this non- Abelian formulation is that we do not need to sum up all Chern 
numbers of the filled one particle states. We just need to evaluated the matrix values field 
strength. This is especially useful for the graphene in a weak magnetic field. Since we are 
mainly interested in the Hall conductance when the fermi energy is near zero. Then one need to 
evaluate all Chern numbers of Landau levels with negative one particle states energies (Landau 
levels of Dirac sea) . The final result is expected to be of the order of unity. It implies most of 
the Chern numbers of the filled Dirac sea are canceled. That is, making one error in evaluating 
each Chern numbers causes a serious problem on the final results. 

Also we mentioned here a numerical technique developed in the lattice gauge theory to 
evaluate lattice topological invariants are quite useful for the evaluation of the £7(1) part of the 
non Abelian Berry connection, which is a two-dimensional analogue of the King Smith- Vanderbilt 
formula for the polarization, this is a lattice version of the one-dimensional line integral [33]. This 
non Abelian method is applied for calculation of the Hall conductance of graphene with realistic 
multi using non orthogonal bases [34J. 

4. Bulk edge correspondence 

Based on the decades of studies for topological states in condensed matter physics, it has been 
widely understood that although the bulk of topologically non trivial states is featureless and 
does not show any fundamental symmetry breaking, near their boundaries and defects as local 
perturbation, there are characteristic local physics governed by the edge states. This is the 
bulk-edge correspondence that can be useful to investigate topological phases both theoretically 
and experimentally [6j. As for the graphene, there have been two types of interesting boundary 
phenomena. One is the boundary states at the zigzag edge with/without magnetic field and 
the other is quantum Hall edge states of the Dirac fermions. Both of them are fundamentally 
controlled by the bulk which we demonstrate here. 

To investigate, let us put the graphene on the cylinder (Fig|6]) with zigzag-bearded edges. 
Let us then take a momentum representation only in tangential e 2 direction as c a (j) = 
I ^ e%k2 ^ c a{ji^2)- Now the total hamiltonian is decomposed as H = J ^H cy i(k2) by a 
one dimensional hamiltonian H cy i(k2) with k 2 as a parameter. H cy i = [t, (ji)cl(ji + l)c (ji) + 

to.{h)cl{ji)c.{j l ) + h.c] where c a (ji) = c a (ji, k 2 ), t. (Jl) = t and t .(ji) = t(l + e^ 2 ^). 
The spectrum of H without magnetic field is shown in FigJTJ There is characteristic zero energy 
states which are localized near zigzag and bearded boundaries [llj. 

One particle Schrodinger equation for each k 2 , H cy i(k 2 )\E(k 2 )) = E(k 2 )\E(k 2 )) by \E) = 

X^i(^»(ii) c »0i' ^2) + ^o(ii)cJ(,7i ; &2))|0) is reformulated using a transfer matrix as ((f) = p/q, 
p and q are mutually prime.) 

v mq£) ) ~ v v>o(o) ; 

where M is a 2 x 2 matrix (see ?? for the precise definition). It implies that the boundary 
condition for zigzag-bearded boundaries is then ip (q£) = (bearded) and V'o(O) = (zigzag). 



Figure 6. Graphene on a cylinder with zigzag 
and bearded edges. Unit cell of the dimensional 
system for the momentum representation by ki is 
shown. 



Figure 7. One particle energy spectrum of the 
graphene without magnetic field labelled by ki- 
There are zero mode boundary states. The red 
ones are localized at the zigzag boundary and the 
blue ones are localized at the bearded one |llll55] . 



Then the spectrum tiifoi) of the edge states is determined by Mixing) = and their position is 
specified as 

for zigzag-bearded: |Mii(q)| > 1 : right j\ oo : (bearded edges) 

< 1: left ji ~ 1 : (zigzag edges) 

When one consider the bearded-zigzag boundaries, the condition is i^ 9 (qi + 1) = (zigzag) and 
= (bearded). Then the spectrum is determined by Mizfa) = and 

for bearded-zigzag: |-M22(eg)| > 1 : right j\ oo : (zigzag edges) 
\M%2 {u) I < 1 : left Ji ~ 1 : (bearded edges) 

4-1. Zero mode boundary states 

Graphene with zigzag edges without magnetic field has a characteristic boundary states 
(Fujita states) which has been also confirmed by a realistic first principle calculation and also 
experimentally [TT) [36j I12| . By the general discussion above, we have {<p = 0) the zero modes 
localized at the zigzag edge when \M\\\ = |1 + e tk2 \ < 1 (2tt/3 < < 4-7r/3) and at the bearded 
edge when \M U \ = |1 + e ik2 \ > 1 (0 < k 2 < 2vr/3, 4vr/3 < k 2 < 2vr)[Hl ES] as shown in FigH 
Actually the existence of the zero localize states is guaranteed by the bulk [35]. Let us first define 
the Berry phase (Zak phase) 7(^2) f° r each k 2 using a bulk hamiltonian h(k) as 

Jz,B(k2) = ~i I dk\ ^r ZB d^ z ,B = or ir (mod 27r) 
Jo 

where ipz,B is an eigen state of the hamiltonian of the periodic systems hz b (compatible 
with the zigzag or bearded edges. This quantization as Z 2 Berry phase is due to the chiral 



symmetry [35, 29] . The zero modes which is compatible for the unit cell of the periodic system 
only exist when 7(^2) = 7T - 

(u \ _ / n ( m °d 27r) : zero modes exist at the zigzag (bearded) edges for &2 
lz,B{ 2) 1 ( mo d 2ir) : no reason to have the zero mode for Jt2 

The condition for the zigzag boundary states, \M\\\ = |1 + e~ lk2 \ < 1, is the condition for the 
Az(k) = 1 + e~ lk2 + e~ tkl encloses the origin by changing k\ :— > 2ir. Since the curve is a unit 
circle centered at 1 + e~ tk2 , the distance between the origin and center |1 + e~ lk2 \ needs to be 
less than unity for the origin to be enclosed. It establishes a bulk-edge correspondence of the zero 
mode boundary states for the zigzag edges. 

With a magnetic field, there still exist zero mode boundary states though the n = Landau 
level also co-exists. Even though the bulk and the edge are both at the same energy, the local 
density of state at E = shows characteristic behavior which is consistent with the localized 
boundary modes [37J. 

4-2. Quantum Hall edge states of Dirac fermions 
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Figure 8. One particle energy spectrum of the 
graphene on a cylinder. The blue (red) lines are edge 
states localized near the left (right) edg es. ky is a 
momentum along the edges and the a xy 's when the 
fermi energy sit there are shown in unit —e 2 /h. The 
numbers are the Chern numbers of the bulk and, at the 
same time, numbers of the blue (red) lines in the gap 
((f> = 1/51). (See also [2]) 



Anomalous quantum Hall effect of graphene as of the Dirac fermions is also discussed with 
the boundaries on the cylindrical geometry. Then counting the light and left edges modes 
between the energy gap where the fermi energy lies, one can assign the Hall conductance using 
the Laughlin argument. Using the transfer matrix formulation given before, we can assign a 
winding number of the edge state on the complex energy surface, which is, in our case, a Riemann 
surface of genus g where g is a number of energy gaps (Landau gaps). Using the setup, one may 
prove a direct relation between the two topological numbers, the one is the Chern number and 
the other is a winding number of the edge state [381 El 05] • It reads physically 

bulk _edge 
u xy %y 

That is, the Hall conductance described by the bulk and the edges are the same. It should be 
like this. This is the bulk-edge correspondence of the Dirac fermions in graphene. In FigJSl 
energy spectrum with the zigzag-bearded boundaries are shown as a function of /C2 and the 
corresponding Hall conductance. It confitms the principle. 
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Appendix A. Landau Level with anisotropic mass 

Let us summarize the standard Landau quantization of electrons with parabolic dispersion 
with (effective) anisotropic mass described by the following hamiltonian H = tv^^-Eltt, 

(rot A = Bz) where 7r = ( ^ x ) , 7Tj = pi — eAi = tt\ and 3l = ( f x ^ y J is a real symmetric 

V % / V *» x v J 

anisotropy matrix. It satisfies [t^x{^b/K), 7r y (£B/K)] = i ( £b = \/~7})- Since the matrix Sl 



is real symmetric, it is diagonalized by the orthogonal matrix as Sl = V^SpV = V^SpV 

S D = diag(6^y), UZy = detH L , & + £ Y = Tr3 L . V=( ^ Y 3 9 G R. Then 

we have H = ltfa D Il, U = Vir, [U x (£ B /h),U Y (eB/h)} = i. 

Now defining a bosonic operator, ([a, a'] = 1) a = {£b/K){JVx + *IIy)/\/2, the hamiltonian 

is written as H = ^f- (f;x{a + o)) 2 — £y( a — a^) 2 ^ where oj = — is a cyclotron frequency. 

Now we define a new bosonic operator ([6,6^] = 1) as a = ub + v*b^ requiring [a, at] = 
[ub + v*tf,u*tf + vb] = \u\ 2 - \v\ 2 = 1. Here noting that £x( a + a" 1 ") 2 - £y(a ~ at ) 2 = 
b 2 {£,x(u + v) 2 - ^y(u - v) 2 } + /i.e. + (btf + + «| 2 + Cy\u - v\ 2 }, we choose as 

Cx(u + v) 2 = £,y{u — v) 2 , u + v = Cy/fr, u — v = ±C\/£x. Assuming £x>£y > an d imposing 
\u\ 2 -\v\ 2 = 1, we have |C| 2 = 1/y 7 ^ = l/(det E L )W, u = i *ggf^, v = -fgffi . 

Finally the hamiltonian is written as 



H = hu(btf + 6t 6)|C| 2 (^y) =/SwH(&t&+i) 

where w= = w\/det Sl = ^v / £x£y- As for the Landau degeneracy, we stressed that it is 
independent of the anisotropy since the magnetic length is independent of the mass. 



Appendix B. Some notations of differential form 

In this appendix, we supplement some details for the differential forms which may not be so 
popular as the standard vector analysis. 

• We have omitted the wedge product (A) which is anti-commuting and useful for the integral 
over the oriented surface. For example, dej) 1 Ad(f> 2 = —d^Adcf) 1 and d^Adcf) 1 = —d^Adcf) 1 = 
0. It implies 

(dG\dG) = {dG\ A \dG) = d^{d^G\ A d(j) u \d v G) = d<f> 1 A d<f> 2 [(diG\d 2 G) - {d 2 G\d l G)] 

• For any (matrix valued) one form X = X^dcf)^, TrX 2 = 

TrX 2 = TrX^X^d^ = TrX^X^d^ = IVX^X^-dftlf ) = -TrX 2 

• As for a matrix valued function u, du^ 1 = —u~ l duu~ 1 . It obeys from = dE = 
d(u)u>~ 1 ) = (duj)u)^ 1 +u>duj~ 1 . 



• Gauge transformation of the non-Abelian gauge strength Fm = A Am + A 2 M where 
Am = u>~ l A' M u + u>~ l du>. 

Fm = d(u} _1 A' m uj) + d(cj _1 dcj) + cj~ 1 A' m 2 cj + u~ l A' M du + u~ l duu~ l A' M u + (cj -1 (iu;) z 
= — u>~ l du>u}~ 1 A' M u) + aj _1 dA' M aj — uj~ l A' M (k>} — ui^ 1 duiui^ 1 du> 

+u~ l A' M 2 u + u~ l A' M du> + u~ l duu~ l A' M u> + (cj -1 ^) 2 = u~ x F' M u> 
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